Macroscopic two-state systems in trapped atomic condensates 
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We consider a macroscopic two-sate system based on persistent current states of a Bose-Einstein 
condensate (BEC) of interacting neutral atoms confined in a ring with a weak Josephson link. We 
demonstrate that macroscopic superpositions of different BEC flows are energetically favorable in 
this system. Moreover, a macroscopic two-state dynamics emerges in the low energy limit. We also 
investigate fundamental limitations due to the noise inherent to the interacting BEC of Josephson- 
ring geometry. We show that the coherent macroscopic dynamics is readily measurable for an 
experimentally accessible range of parameters. 

PACS numbers: 03.75.Kk, 37.10.Gh, 85.25.Cp 



Realization of macroscopic quantum two-state systems 
has been a challenge for cold-atom BEC physics from 
early stages of its experimental development. The re- 
search in this area has been primarily focused on many- 
particle dynamics within double-well trapping potential 
At the same time, despite an appealing similar- 
ity with the microscopic single particle two-state system, 
this geometry provides no easy way to achieve a superpo- 
sition of distinct many-body quantum states. Indeed, a 
collection of non-interacting (or weakly interacting) bo- 
son atoms confined in a double-well trap condenses into 
a BEC with the "product" wave function not suitable 
to form a macroscopic two-sate configuration. Repul- 
sively interacting particles in such trap favor [s| Fock 
states — the system enters Mott or "Coulomb blockade" - 
like regime. Hence, attractively interacting particles be- 
come the only option in such geometry. In the latter case, 
atoms correlate jjj forming a "Schrodinger cat" state 
*x<o = IIL M^) + IL M^i)]/V2, where ^ L/R (r) dc- 
note the single particle states localized the left /right well 
of the trapping potential. Observation of this macro- 
scopic superposition, however, is extremely challenging. 
The only manifestation of a coherent superposition is 
the presence of the off-diagonal matrix elements in the 
macroscopic two-state basis. In the case of the double- 
well trapping potential these matrix elements are pro- 
portional to the probability of all N particles to tunnel 
through the barrier, which is extremely small. 

In this paper we study a cold atom based macroscopic 
two-state system (a qubit) based on a persistent current 
BEC- Josephson system [5j. The two macroscopic states 
are metastable current-carrying states of a BEC confined 
in a Josephson ring trap [6|, |7|. Such systems have be- 
come experimentally available due to recent successes in 
dynamical BEC trapping Q. We start by deriving an 
effective Schrodinger equation describing the low-lying 



energy states of the system in terms of the phase dif- 
ference across the Josephson junction. For that we in- 
troduce a simple anzats that parameterizes the low-lying 
states as a superposition of the condensate states with 
different phases. Then we consider a better ground state 
anzats that accounts for the fluctuations in the systems 
arising from the interparticle interactions. We briefly 
discuss their influence on system's dynamics and then 
analyze their effect on the detection of the systems's 
state. Particularly, we show that the admixture of "non- 
condensate" particles acts as an effective noise in the 
time-of-fiight (TOF) images of the BEC and therefore 
imposes limitations on the resolution of such measure- 
ments. We briefly discuss these fundamental limitations 
and argue that weakly coupled Bose systems provide a 
good candidate to observe macroscopic quantum tunnel- 
ing (MQT) and macroscopic quantum superpositions. 

We consider a system of N locally interacting bosons 
confined to the rotating external potential V correspond- 
ing to a ring with a thin cut (a barrier) . We assume that 
the transverse dimension of the ring is small compared 
to the healing length of the bosons and therefore the sys- 
tem can be considered effectively one-dimensional. The 
Hamiltonian of such system can be written as Q 
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where rotation frequency of the ring (i.e., the barrier), 
R is the radius of the ring, A = Aira/Sm, where a is 
scattering length and S is the ring's cross section area. 
Here and in the following we will use units with H = 1. 

We start by evaluating the energy of the ground state 
of the Hamiltoian ([IJ using Gross-Pitaevskii approach. 
That is, we assume that the ground state wavefunction 
is a product, ^(ri, ...,r/v) = x( r i)---x( r N) and minimize 
the the functional E = J dri...dr^'^*H' i k by varying 
it with respect to the single-particle state \- It is ob- 
vious, however, that, since the system is homogenous 
everywhere except in the small region at the barrier, 
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const or x ~ e «*(r) (except in the vicinity 



of the barrier). Moreover, the phase-dependent terms 
in E are ~ / dr[(V$) 2 /(2m) + fiEV$], and therefore $ 
must be a linear function of distance along the ring, i. 
e., x ~ e l ^ 6 '/ 27r , where is the phase difference across 
the barrier and 9 is the azimuthal angle parameterizing 
position along the ring. Evaluation of E for such prod- 
uct state yields Ncj> 2 /2mL 2 - NCI^/2tt + XN(N-1)/2L, 
where L is the circumference of the ring. This expression, 
obviously, does not account for the contribution due to 
the barrier, i.e., the Josephson energy. At the barrier the 
particle density, N\x(r)\ 2 is strongly dependent on r and 
therefore must be calculated self-consistently. It can be 
shown, however, that to a good approximation, contribu- 
tion of the barrier region into the system's energy can be 
cast in the form — Ej cos A, where the Josephson energy 
Ej is independent on |9|. Thus we find that 

E((j>) = XN(N-l) /2L +N(<j) -<j) ) 2 /2mL 2 -Ej cos (f>, (2) 

where </>o = mL 2 il/2ir. 

When (f>Q = it the effective potential is a symmetric 
double well, corresponding to two macroscopically differ- 
ent states, i. e., carrying different persistent currents. 
One can see, however, that such degeneracy is lifted 
by quantum fluctuations, leading to macroscopic quan- 
tum tunneling. In order to see this, let us evaluate the 
system's energy for the superposition state C J d<j)^^\ 

where is the above product state and C is normaliza- 
tion constant. A straightforward calculation shows that 
the expectation value of the interaction energy, i.e., of 
the last term in Eq. ([1]), is lower by XN/2L than that for 
the localized (in <fi) state ■ 

We note, however, that while the interaction part of 
the Hamiltonian in Eq. (p} favors superposition, the first 
two terms , obviously, "prefer" the localized state. There- 
fore we search an optimal ground state wavefunction in 
the following form [1CJ, 1 : 



#(!•!,..., rjv) = v^V/24^ / d^)¥i 0) (ri,...,rjv), (3) 



where tp(<p) i s to be defined by minimization. In or- 
der to evaluate expectation value of the energy for the 
wavefunction in Eq. ^ we note that the states 
are approximately orthogonal. Indeed, 5n{4> ~ 4>') = 



N 



J dr 1 ...dr N ^>*^/ yields S N (</>) = [sin(0/2)/0/2] 
e -~N<f> /24 This is a rapidly varying function — the combi- 
nation J N /2Att8n{4>) approaches the true 5- function for 
N 3> 1. Moreover, one can see that for any few-particle 
operator A the following identity holds: 



<*S )) |i|*w> = ^W-^)<*7|i|*J k UJ )[i+0(i/JV)]. 

As the result, one could think that the expectation value 
of ([!]) calculated with respect to the wavefunction §S§ 
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is merely / #|^)| 2 ^ 0) |#|*f) = / d^(<f>)\ 2 E(cb). 
This, however, is not the case. The expectation value 
of the interaction part of the Hamiltonian is amplified by 
iV 2 and the neglected 0(1 /N) terms become important — 
the terms appearing due to the finite width ol Sn(4>) con- 
tribute to the first order in N. A more careful evaluation 
of the interaction expectation value gives 



XN(N— I) 
2L~ 



##V(</>0'W)' 



/247T iVW W tan^ V ; 



The finite width of 5n{4>' ~ 4>) is clearly non- negligible. 
We use the identity N 2 (j) 2 8 N {(j>) / 12 = 12c§5 N (<t>) + 
N5n(4>) an d evaluate the integral over ip containing 
dtpdtpSjy ((j)' — 4>) by parts, finally arriving at the result 
correct to @(N) 
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XN 
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-E(4>) 
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The first term in the right hand side of Eq. ^ is the 
negative offset mentioned earlier. The second term is a 
positive gain in energy due to variation of <j). Together 
with the last term in Eq. ([5]) it can be viewed as kinetic 
and potential energies of a "phase-particle", whose dy- 
namics obeys Schrodinger equation Hegi]) = Etp, where 
H cS = -(QX/L)d 2 +E(cp). 

The ground state for such "particle" is the symmetric 
superposition of two states in each well of the effective po- 
tential E ((/>), separated from the next excited state, i.e., 
the antisymmetric combination, by the tunnel splitting 
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exp [-r]^S(j) 3 /72} 



4>r~4>l, (7) 



where the exponential prefactor is the frequency of small 
oscillation in each well, 4>l and <pR are positions of 
each minima, rj — \J NS/inaL is Tonks parameter 12 1 
(typically 77 = 10 — 100 for cold atom systems), and 
a = mEjL 2 /N . When potential barrier is small, a ~ 1 
[lH . Energy AE defines the timescale of tunneling tran- 
sitions between states with two different persistent cur- 
rents. 

The states with different persistent currents can be 
readily detected by looking at absorbtion images of the 
density distributions in time-of-flight (TOF) measure- 
ments @,Hj], e.g. Fig. (1). For long TOF f single par 

X4> 



ticle wavefunctions yT (r) become Fourier transforms 



of the initial single particle wavefunctions, x$ ( r ) ~ 
X</>(Q)lq=2mr/t ■ The outcome of the measurement is best 
understood by considering different moments of p(q). A 
straightforward calculations yields 



(p(q)) = J ^|^)| 2 P0 (q) (8) 
(p(q)p(q')) = / d<t>\^)\ 2 PM)PM) 
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FIG. 1: (Color online) Expected TOF images of the 
Josephson-ring system in a double-well regime, (a) TOF den- 
sity at the center as a function of (j>. Typical double-well 
potential minima 4>l, 4>r are shown to give relevant scale on 
<j>. The typical TOF image corresponding to measurement 
outcome (b) <j>L and (c) 4>r. In both cases the cross shows 
the center of the image. 



where ^(q) = (N / S)\x<f>{<i)\ 2 ■ These moments corre- 
spond to a stochastic process: the probability to observe 
the outcome p0(2mr/ixOF) and, hence, a particular value 
of </>, is given by |^/>(<^>)| 2 . Therefore, the measurement 
"chooses" a single term in the sum ©, and the function 
ip(4>) becomes a quantum- mechanical wave function of a 
macroscopic object — the superfluid current of the entire 
condensate. 

The above picture obviously is an approximation: We 
have assumed that all particles are "in the condensate" , 
e.g., Eq. ([3]). For weakly interacting bosons the ground 
state is no longer a product, but the well-known Bogol- 
ubov state. Therefore we consider a modified ground 
state wavefunction, still in the form of Eq. @ , but with 
modified i.e., 



*i 0) ) 1*0) ~ ( 



Et t 
C n a <f>+2iTn a <l>-2nn< 



)" /2 |0>, (9) 



where a 



(t) 



is a creation (annihilation) operator corre- 



sponding to a single-particle state X<Lui r ) = X<j>( r ) ein6 '> 
and c„ are variational coefficients [14]. By choosing 
the trial wavefunction in the form of Eq. ^ we as- 
sume that the dynamics is adiabatic: The slow collec- 
tive variable <f) is coupled to fast degrees of freedom, i.e., 
the fluctuations. Indeed, due to the finite size of the 
the system the quasiparticle spectrum is discrete with 
gaps ~ (speed of sound/i) = N/(mL 2 r]), which is much 
greater than the energy scale associated with the tun- 
neling between two wells, e.g., Eq. J7|. As a result the 
quasiparticles remain in the ground state and readjust to 



the variations in 0. 

Corrections to the effective energy functional, e.g., 
Eq. ©, turn out to be small in the limit of weak inter- 
actions. The new ground state energy can be evaluated 
by using Eq. which, as can be directly verified, still 
holds for the modified wavefunction, e.g., Eq. ©. Then, 
evaluating (^rplH^^} and varying c„'s (see Ref. 0] for 
similar calculation) we find that corrections to the func- 
tional E in Eq. (J6j) arc suppressed b y fa ctor 1/r), .i.e., are 
small in the weak interaction limit [151 ] . 

The effect of fluctuations, however, turns out to be 
quite appreciable as far as measurement of close current- 
carrying states is concerned. Indeed, while the number 
of non-condensate particles is small, the difference be- 
tween density profiles (in the TOF measurements) corre- 
sponding to different values of cj>, i.e., different persistent 
currents, is also small. In the opposite case the barrier 
is to high to allow for tunneling on a reasonable time 
scale. Moreover, as we will see, both the tunneling ex- 
ponent in Eq. (0 and the signal-to-noise ratio due to 
" non-condensate" particles (to be defined below) are con- 
trolled by the same parameters. To see this let us eval- 
uate the set of correlation functions as in Eq. but 
for the modified wavefunction given by Eq. ^ . In doing 
so we again use Eq. (0}, and so the calculation reduces 
to the evaluation of averages (fy^lp^p^')...] 1 ^^), where 
TOF density operator p(q) = Y^ n ,n> X*<j,,„(<i)X<t>,n' a„' 
(again X<t>,n{o) is the Fourier transform of X<f>,n(f))- The 
first two moments arc 



(P(q)> 
(p(q)p(q')) 



d^(4>)\ 2 ~ P ^) 

#|V#)| 2 ^(q)^(q') 



(10) 
A (q,q')] 



where p^q) = P<j} (q) + En#o( a 0,n a <M>IX0,n(q)| 2 is the 
renormalized density and 

A^(q,q') = ATX^,n(q)X^n(qOH+2,rn a <H-27ro) (H) 
n 

(q)x<^n(q')( a 0+ + c. a . 

In evaluating Eqs. (fTUl ITT]) we have used Wick's theo- 
rem as well as the fact that in the limit of weak inter- 
actions (al +27Tn a^ +2ml )\ n=0 ~ N > (4, +27m a 0+2lm )|„^ o 

and therefore products {a^ +2vr a^ +2 ^n){a\ +2v n' a <l>+^r l ') , 
etc., with n, n! ^ 0, can be neglected. Moreover, since 
A^(q, q') is small compared to p^(q)p^(ci'), we can re- 
place it by A ff (q, q') in Eq. (flTJ|) . - we are interested in 
the situation when the barrier separating the two per- 
sistent currents is small, 8(f> <C 7r. As a result A 7r (q, q') 
can be viewed as the correlation function of noise super- 
imposed with "signal" ^(q). It is natural, therefore, to 
introduce a signal to noise ratio (SNR) as 



SNR 



[Pn-6<t>/2(0) ~ Pn+S(f,/2(0)Y 

A w (0,0) 



(12) 
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where we define the strength of the signal as the differ- 
ence in the TOF particle densities at the center for two 
realizations with phases %±5tfi, e.g., Fig. 1. The SNR in 
Eq. (fT2|) can be easily evaluated: At the center X0,„(q = 
0) = B/((j)/2% + n), where B = (1/2%) J dr X ^(r).' Then 
we obtain that 



\ a ti a n) , (a n a- n ) . 



l-4n 2 (l+2n) 



(13) 



where a shortcut a! = a S+27rn nas been used to re- 



duce notations. The averages in Eqs. (|T3|) can be 
found from the normal and anomalous Green's func- 
tions of the interacting boson system, i.e. (a' n a n ) = 
f t£G(ti, ioj) and (a„o_„) = J" |^.F(n, iw), where 

G(n,iu) = aild Jr ( n '* W ) = u' X +(et)' ' With 

e% = h 2 ((j)- cj) Q + 2%n)/2mL 2 + Xp/2 and = 

y (en) 2 — (^y°/2) 2 [HI ■ After straightforward calculation 
we find that 



SNR ~ 7.84ri<5</> 2 



(14) 



Eq. (|T4f is the principal result of this paper. It re- 
flects the fundamental difference between the phase dy- 
namics of the Josephson qubits and single body quan- 
tum mechanics: While quantum mechanics, in princi- 
ple, allows one to measure particle's coordinate with in- 
finite precision, the accuracy of the phase measurements 
in Josephson qubits is limited by their many-body na- 
ture. Remarkably, the precision for such measurement is 
controlled by a single parameter 77, i.e., the interaction 
strength. Moreover, comparing Eq. ([7]) with Eq. (| 14[) we 
see that both the tunneling exponent and the SNR are 
controlled by the same parameters. This fact is not sur- 
prising: while stronger interactions enhance MQT, e.g., 
Eq. 0, they decrease the SNR due to the suppression 
of the condensate density for stronger scattering. The 
tunneling exponent, however, increases faster (in abso- 
lute value) with the growth of 8<j) and therefore to keep 
it small together with the condition SNR > 1 we need 
77 3> 1. Thus coherent MQT is observable only in weakly 
interacting BEC systems. For typical BEC experiment 
77 - 10-100 0," 
in (0 equals to 



17] . If we keep the tunneling exponent 
1 we obtain SNR - 300 - 600. 
To conclude, we have demonstrated that a cold atom 
Josephson-ring system can be described by an effective 
single-particle Schrodinger equation with an effective po- 
tential that can be controlled by the system's rotation. 
We analyzed the statistics of the TOF measurements that 
allow one to determine the system's current, i.e., phase 
across the Josephson junction. We found that fidelity of 
such measurements is limited by the ground state fluctu- 
ations and is controlled by the Tonks parameter. 

Finally we conjecture that similar limitations are likely 
to take place in superconducting flux qubits. While nu- 



merous studies of such systems have been carried out, in- 
cluding the derivation of the effective Schrodinger equa- 
tion, etc., generalization of the results obtained in this 
paper to the superconducting case, however, does not 
seem to be straightforward due to a different nature of 
the ground state for a fermionic system. Thus we believe 
that the extension of our results to superconductors is an 
interesting direction for future work. 
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